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We argue that lattice calculations of the r[ mass in QCD with N c = 2 colors performed at non¬ 
zero baryon chemical potential can be used to study the mechanism responsible for the mass of the 
rf. QCD with two colors is an ideal laboratory because it exhibits confinement, chiral symmetry 
breaking and a would-be U(1)a Goldstone boson at all densities. Since the instanton density and 
the confinement scale vary with density in a very different way, instantons are clearly distinguishable 
from other possible mechanisms. There is an instanton prediction for the rf mass at large density 
that can be compared to lattice results. The density dependence of the instanton contribution is a 
simple consequence of the integer topological charge carried by the instanton. We also argue that 
N c = 3 color QCD at finite isospin density can be used in order to study the origin of OZI-violation 
in the scalar sector. 


I. INTRODUCTION 

The U{1)a puzzle in QCD is related to the absence of 
a ninth Goldstone boson connected to the spontaneous 
breakdown of the 17(1 )a chiral symmetry jlj]. It was re¬ 
alized soon after the discovery of QCD that the U(1)a 
symmetry of the QCD lagrangian is anomalous, but it 
was also noted that the divergence of the U(1)a current 
is itself a total divergence j|. Superficially, it would then 
seem that the 17(1 )a anomaly is not sufficient to remove 
the U(1)a Goldstone boson. The puzzle was resolved 
after Belavin, Polyakov, Schwartz and Tyupkin discov¬ 
ered topological structures, instantons, in QCD [|J. ’t 
Hooft showed that instantons lead to the violation of 
axial charge Q, and that instantons induce an effective 
(27Vy)-fermion operator which contributes to the rf me¬ 
son mass 

Since then, lattice QCD calculations have identified 
instantons, verified the presence of fermion zero modes, 
and established their relation to the mass of the rf |ij. 
Also, the instanton liquid model was expanded into a 
phenomenologically successful description of chiral sym¬ 
metry breaking and the 17(1) a anomaly 0 |. Unfor¬ 
tunately, we cannot compute the instanton contribution 
to the rf mass from first principles. As a consequence, 
there are still speculations that the rf mass is related to 
structures with fractional topological charge that do not 
appear in the classical limit, or that the rf mass is in 
some way related to confinement ||. The latter sugges¬ 
tion was first made by Kogut and Susskind prior to the 
discovery of instantons in QCD ||lof . 

Even if we cannot do a parameter-free calculation of 
the rf mass in QCD we can still try to distinguish dif¬ 
ferent mechanisms for generating the rf mass by their 
scaling behavior. In QCD, of course, any contribution 
to the rf mass has to be proportional to the QCD scale 
parameter A qcd- Witten suggested that the number 
of colors N c could be used as a parameter jn|. He ar¬ 
gued that the instanton contribution scales as exp(— N c ) 
whereas effects related to confinement give ?tt 2 , ~ 1 /N c . 
However, the relation to 2 , ~ exp (—A?) is only correct for 


very small instantons, p -C A qcd- Indeed, we recently 
argued that the instanton contribution to the rf mass 
also scales as ?n 2 , ~ l/iV c Q. 

Another possibility is to use the temperature T as a 
parameter. The density of instantons is expected to be 
suppressed by a large power of T, (N/V) ~ T b ~ 4 , where 
b = llN c /3 — 2Nfj3 is the first coefficient of the QCD 
beta function Effects related to confine¬ 

ment or topological objects other than instantons would 
presumably have a different dependence on temperature. 
The problem with this idea is that the suppression of in¬ 
stanton effects is related to perturbative color screening. 
This implies that the power law suppression only applies 
if the temperature is larger than the critical temperature 
T c for chiral symmetry restoration jlTj]. In order to have 
rigorous theoretical control over instanton effects we have 
to consider the quark gluon plasma phase instead of the 
hadronic phase. 

In this note we suggest using the baryon chemical po¬ 
tential p as a parameter. We shall argue that this pa¬ 
rameter is more useful than the number of colors or the 
temperature because there is chiral symmetry breaking 
and a hadronic phase for all values of p, and there is the¬ 
oretical control over both the instanton contribution to 
the rf mass and the scale of confinement effects. We first 
formulate our proposal in terms of SU( 2) gauge theory 
at finite baryon density. We then show that a similar 
situation arises in SU( 3) gauge theory at finite isospin 
density. Both of these theories can be studied with lat¬ 
tice algorithms that are available today. 

We should note that many of the features that we shall 
discuss also apply, with some modifications, to SU( 3) 
gauge theory at finite baryon density JlJ]. However, this 
theory cannot be studied on the lattice at present, and 
we shall not discuss it in detail. 


II. QCD WITH TWO COLORS 

Let us summarize some of the salient features of SU (2) 
gauge theory at zero and non-zero baryon chemical po¬ 
tential [|l9|, |o], fn[ HJ. For simplicity, we will con- 
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FIG. 1: Phase structure of N c = 2 QCD at finite baryon 
chemical potential and temperature. We consider the diquark 
phase m = 0,j —> 0. On the T yf 0 line the symmetry break¬ 
ing pattern is 5(7(4) —+ 5p(4). Above T c , the 5(7(4) symme¬ 
try is restored. At fj, ^ 0 the 5(7(4) symmetry is explicitly 
broken to 5(7(2) x 5(7(2) x (7(1). For T < T c this symmetry 
is spontaneously broken to 5(7(2) x 5(7(2). The dashed line is 
the location of the deconfining phase transition. There is no 
strict order parameter for this transition at low density and 
the transition may just be a rapid cross-over. 


centrate on Nf =2 flavors. 5(7(2) gauge theory has a 
meson spectrum which is very similar to 5(7(3) QCD. 
Baryons, on the other hand, are bosons rather than 
fermions and their spectrum is very different as com¬ 
pared to N c = 3 QCD. Because the 5(7(2) gauge group 
is pseudo-real, there is a Pauli-Gursey symmetry which 
relates quarks and anti-quarks. This symmetry mixes 
the quark-anti-quark condensate ( q a q a ) with the diquark 
condensate (e ab q a T C^$T 2 q b ). Here, a, b are color indices 
and t 2 is the anti-symmetric 5(7(2)* flavor matrix. As 
a result the chiral symmetry group is 5(7(4) rather than 
SU(2) l x 5(7(2)*. 

Let us consider chiral symmetry breaking in the pres¬ 
ence of a small source term 


anti-quark condensed phase. If the baryon chemical po¬ 
tential is non-zero the 5(7(4) symmetry is broken ex¬ 
plicitly. In the following, we shall consider the diquark 
phase m = 0, j —> 0. In this case, the phase diagram 
is particularly simple, see Fig. 1. Chiral symmetry is 
spontaneously broken at T = n = 0. There is a critical 
temperature T c such that chiral symmetry is restored for 
T > T c . However, most likely, there is no phase transi¬ 
tion as a function of the chemical potential for T < T c . 
This, of course, is the main feature that distinguishes the 
5(7(2) theory at non-zero /i as a laboratory for studying 
the (7(1 )a anomaly. 

The effective lagrangian for the singlet pseudoscalar 
Goldstone boson is 


£ = fp [( d 0 (j)) 2 - v 2 (d,4>) 2 ] - V{(j>). (2) 

The decay constant and Goldstone boson velocity can be 
determined in perturbation theory. At leading order, the 
result is [|d| 


fp = 




(3) 


The potential V (</>) receives contributions from instan- 
tons. We find V (0) = —Ap cos {(j)+6) where 9 is the QCD 
theta angle. If the chemical potential is big, p, A qcd, 
large instantons are suppressed and the coefficient A can 
be determined in perturbation theory. We find |^4|, ^5| 


Ap = C 2t2 6n^ 





A" 2 (4) 


with 


_ 0.466 exp(—1.679 N c )1.34 N t 
CN °' Nf ~ {N c — 1)!(-ZV C — 2)! • (5) 

At large ^ the gap A can also be determined in pertur¬ 
bation theory. We get |26|. |28|, ||i| 


C s = mq a q a + - ( e ab q aT C'y 5 T 2 q b + h.c.) . (1) 

Chiral symmetry is broken according to 5(7(4) —> 5p(4) 
in the low-temperature and low-density phase for both 
j = 0,m —> 0 and m = 0, j —> 0. In the case m ^ 0 
the order parameter is the quark-anti-quark condensate 
( q a q a ). There are 5 Goldstone bosons, three pions n, the 
scalar diquark 5 and the scalar anti-diquark 5. Because 
of the U(1)a anomaly the would-be singlet Goldstone bo¬ 
son, the r]' , is heavy. In the case j / 0 the order param¬ 
eter is the diquark condensate (q T C^T 2 q). Again there 
are 5 Goldstone bosons, three pions tt, the sigma a and 
the scalar diquark 5. The would-be singlet Goldstone 
boson, the pseudoscalar diquark P, is heavy. 

For n = 0 all directions of the source term (rn. j) = 
TOo(cos(a),sin(a)) are physically equivalent. In partic¬ 
ular, the masses of the pseudoscalar diquark P in the 
diquark phase is equal to the mass of the r( in quark- 


A = 


5127r 4 f>Q/n7 5 




( 6 ) 


where the parameter b' 0 = exp(— (tt 2 +A)(N C — 1)/16) con¬ 
trols the size of non-Fermi liquid effects j3(], |l|. Using 
equ. (||) we can determine the mass of the pseudoscalar 
Goldstone boson. We have 



(7) 


The result has the structure of the Witten-Veneziano re¬ 
lation where Ap plays the role of the topological suscep¬ 
tibility. We note, however, that xtop = 0 as expected 
for a theory with massless fermions. In fact, Ap governs 
local, not global, fluctuations of the topological charge. 
For a dilute gas of instantons A P = (N/V) where (N/V) 
is the density of instantons. This relation is exact in the 
limit of large baryon density. We also note that in the 
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law is directly related to the topological charge of the 
instanton. Contributions from instantons with charge 
two or larger are suppressed by additional powers of 
{Aqcd / rf)- By the same token we expect the contri¬ 
bution of hypothetical objects with fractional charge to 
dominate over instanton effects at large /x. 

In practice, lattice calculations have to be carried out 
at finite diquark source or non-zero quark mass (or both). 
If the baryon density is small the pseudoscalar diquark 
and rf meson are very heavy and the effect of a non¬ 
zero quark mass is small. At large g, however, the quark 
mass contribution is more important. If /x 3 > A qcd the 
quark mass contribution to the effective potential can 
be computed in perturbation theory. Using the methods 
described in |32) we find 

V m (4>) = det(M)e-* + h.c., (9) 



FIG. 2: Fig. a) shows the pseudoscalar Goldstone boson mass 
in A c = 2 QCD. The solid line marked by squares is the result 
of the instanton calculation. We have used the one-loop run¬ 
ning coupling constant with A qcd = 170 MeV. The solid line 
marked by circles shows an estimate based on the assumption 
Tripfp ~ cr 2 . The overall scale of this curve is arbitrary. For 
comparison, we also show the energy gap 2A and the con¬ 
finement scale. Again, the overall scale is arbitrary. Fig. b) 
shows the scalar Goldstone boson mass in N c = 3 QCD at 
finite isospin density. 


where M is the mass matrix. We observe that the topo¬ 
logical susceptibility does not vanish if the quark mass 
is non-zero, \top ~ det(M)A 2 . The contribution to the 
mass of the would-be U(1)a Goldstone boson is 


m 


2 _ 
p ~ 


32 A 2 


-m u m d - 


( 10 ) 


In the limit fi —> oo the m q rf 0 contribution will eventu¬ 
ally dominate over the instanton contribution. However, 
even for quark masses as large as m u = m d = 40 MeV 
the instanton contribution is expected to dominate for 
chemical potentials that can be achieved on the lattice. 

We now comment on possible effects related to con¬ 
finement (H], |3^, Q. At large baryon density the gap in 
the fermion spectrum is much larger than the QCD scale 
parameter, A A qcd- Since the diquark condensate 
is a color singlet, there is neither screening nor a Higgs 
effect operating at scales below A. As a consequence, we 
expect the color SU (2) to be confined at all densities. If 
there are effects related to confinement that contribute 
to the mass of the rf, then these effects should persist at 
all densities. 

We do not know how to compute confinement related 
contributions to the mass of the rf. However, if these 
effects exist then they should be governed by the SU (2) 
confinement scale. Rischke et al. observed that the pure 
glue theory below A is characterized by a non-trivial elec¬ 
tric polarizability and magnetic permeability iflB) . The 
effective action is 


instanton liquid model Ap ~ ( N/V ) is very well satisfied 
even at zero baryon density. 

We observe that equ. 0 ) implies 



As expected, the instanton contribution to the U(1)a 
mass is suppressed at large baryon density. The power 


Seff = 27 i d4x ( £jg ° ' ^ " \ §a '^“) ’ (11) 

where we have suppressed higher order terms that are 
suppressed b y Op / A). The main feature of the effective 
action equ. (|ll|) is that the electric polarizability e ~ 
1 + g 2 rf 2 /(18 tUA 2 ) is very large in the limit g —> oo. As 
a consequence, the SU( 2) pure gauge theory described 
by equ. & is confined, but the confinement scale is 
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exponentially small ||i| |3(| 

. . ( 8tt 2 3-6^ , 

A„„,~Aexp(-^ —(12) 

If the mass of the would-be U(1)a Goldstone boson is 
dominated by confinement effects then we expect that 
m 2 P fl ~ cr 2 ~ A 4 on/ . In Fig. 2 we compare this scal¬ 
ing relation with the instanton prediction for the pseu¬ 
doscalar Goldstone boson mass. Of course we do not 
know the constant of proportionality relating mpfp and 
the string tension a. In Fig. 2 we have arbitrarily scaled 
mp ~ A 2 con j/fp to match the instanton prediction at 
fi ~ 500 MeV. We find that the instanton and confine¬ 
ment related contributions clearly scale differently. We 
also observe that the instanton contribution will always 
dominate at large p. 

We should note that the effective theory described by 
equ. (0 is not covariant. As a consequence, the time¬ 
like and space-like string tensions are not the same. We 
can restore covariance by rescaling the time coordinate 
in equ. © according to x' 0 = Xq/ yfe. This implies that 
the two string tensions differ by a factor yfe. We have not 
tried to take this effect into account, because there is no 
definite theory that would predict the dependence of mp 
on the spacelike and timelike string tensions. We should 
also note that the estimate equ. ([ 12 ]) of the confinement 
scale only has exponential accuracy, so that we cannot 
reliably predict possible factors yfe in the pre-exponent. 


At large isospin density there is a U(1)a would-be 
Goldstone boson, the aj. The effective lagrangian for 
the a(J~ is identical to the effective lagrangian for the r( 
in N c = 2 QCD, see equ. (]|). If jij A qcd the pa¬ 
rameters in the effective lagrangian can be computed in 
perturbation theory. The decay constant and Goldstone 
boson velocity are given by 


fs = 



(13) 


The instanton contribution to the effective potential is 
given by Ag cos {<j> + 9) with 


Ag = 6)3.2 


167T 4 r(f) 

235/6 


<16 




35/3 

A’ 2 , (14) 


where the gap M and superfluid density <f>g are given by 
M = 512n%^g~ 5 exp , (15) 

(16) 

The mass of the aj satisfies a Witten-Veneziano type 
relation 


m 


2 

OLQ 


Ag 

2 / 1 ' 


(17) 


III. QCD AT FINITE ISOSPIN DENSITY 


Again, Ag is related to the density of instantons. At 
large /// the oj mass scales as 


The large mass of the rf implies that violations of the 
OZI rule in the pseudoscalar meson sector are substan¬ 
tial. However, the rf sector is not the only channel in 
which OZI violation is large. In particular, the OZI vio¬ 
lating mass difference between the isovector-scalar a 0 and 
isoscalar-scalar a meson is almost as large as the rf — n 
splitting. We have argued that the ao~a splitting is also 
dominated by instantons 0H3- In this section we show 
that this idea can be checked in N c = 3 QCD at finite 
isospin density. 

In QCD with N c = 3 colors and Nf = 2 flavors chiral 
symmetry is broken according to SU(2)p x SU(2)p. If 
the quark mass is non-zero there are three almost mass¬ 
less pions and a heavy would-be U(1)a Goldstone boson, 
the rf . There is a relatively light scalar-isoscalar meson, 
the <7, and a heavy scalar-isovector meson, the a 0 , which 
is close in mass to the 77 '. 

The effect of a non-zero isospin chemical potential term 
gi(u^u — dJd) was studied in Q. If gi > 0 the isospin 
chemical potential favors up quarks over down quarks. 
As a result, the mass of the positive pion is reduced. If 
gi > 771^/2 pion condensation takes place and the chiral 
order parameter starts to rotate from the (qq) direction 
to the (qi"f 5 T~q) direction. There is one exact Goldstone 
boson, the 7r + , and two heavy pions. 



The difference in the power law suppression as compared 
to equ. (||) is related to the difference between the beta 
functions for N c = 2, 3. 

Finally, we observe that the pion condensate is a color 
singlet and the color SU (3) is expected to be confined at 
large isospin density. The confinement scale is given by 

A„„~Aexp(—(19) 

In Fig. 2b we plot the instanton prediction for the mass 
of the ag" at large isospin density. Again, we compare 
the instanton prediction to the scaling relation m ao ~ 
cr/ fs- We observe that the instanton contribution is more 
strongly suppressed as compared to the N c = 2 result. 
As noted below equ. (^ 8 [) this is related to the fact that 
the beta function is larger. The scale a/fs, on the other 
hand, drops off more slowly as compared to the N c = 
2 case. This is a consequence of the fact that the gap 
is larger and as the result the polarizability is smaller. 
As asymptotically large isospin density, however, the a 0 
mass is again dominated by instantons. 
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IV. SUMMARY 

We have argued that SU (2) QCD at finite baryon den¬ 
sity and SU (3) QCD at finite isospin density can be used 
in order to study the mechanism for generating the mass 
of the would-be U(1)a Goldstone boson in QCD. The 
main point is that there is a hadronic phase at all den¬ 
sities, and that both the instanton contribution and the 
confinement scale are calculable at large chemical poten¬ 
tial. Furthermore, both SU{ 2) QCD at non-zero baryon 
chemical potential and SU{2>) QCD at non-zero isospin 
density can be studied on the lattice with presently avail¬ 
able methods [icj . 

We would like to mention some possible difficulties 
with our proposal. We have assumed that there is no 
phase transition along the finite baryon chemical poten¬ 
tial axis in the phase diagram. This assumption is based 
on the observation that the phenomenologically estab¬ 
lished symmetries of the low density phase agree with the 
calculated symmetries of the high density phase. This 
means that a phase transition is not required, but of 
course a transition is not forbidden either. The question 
of whether or not there is a phase transition can be stud¬ 
ied using lattice simulations. We should note that even if 


there is transition that separates the low and high-density 
phase, we can still study the question whether the mech¬ 
anism for generating the mass of the U(1)a Goldstone 
boson is the same in both phases. 

We have used perturbation theory in order to com¬ 
pute the gap, the instanton density, and the confine¬ 
ment scale at non-zero chemical potential. We do not 
know how reliable leading order perturbation theory is 
for baryon chemical potentials that can be achieved in 
lattice calculations. Indeed, it was argued that the per¬ 
turbative expansion for the gap converges very slowly 
and that the instanton contribution is very sensitive 
to the value of the QCD scale parameter J25|. However, 
the power law behavior of the C7(1)^4 Goldstone boson 
mass is quite robust and a simple reflection of the topo¬ 
logical charge carried by the instanton. Furthermore, the 
question whether the U[1)a Goldstone boson mass scales 
with the instanton density, m 2 p fp ~ (N/V), or the string 
tension, m 2 p fp ~ a 2 , can be answered directly from the 
lattice data, without resort to perturbation theory. 
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